Bounding the Resource Availability of 
Partially Ordered Events with Constant Resource 

Impact 


Jeremy Frank 


Computational Sciences Division 
NASA Ames Research Center, MS 269-3 
frank@email.arc.nasa.gov 
Moffett Field, CA 94035 


Abstract We compare existing techniques to bound the resource availability of 
partially ordered events. We first show that contrary to intuition, two existing 
techniques, one due to Laborie and one due to Muscettola, are not strictly compa- 
rable in terms of the size of the search trees generated under chronological search 
with a fixed heuristic. We describe a generalization of these techniques called the 
Flow Balance Constraint to tightly bound the amount of available resource for a 
set of partially ordered events with piecewise constant resource impact We prove 
that the new' technique generates smaller proof trees under chronological search 
wfith a fixed heuristic, at little increase in computational expense. We then show 
how to construct tighter resource bounds but at increased computational cost. 


1 Introduction 

Scheduling is about simultaneously satisfying temporal constraints and resource avail- 
ability constraints. Following Muscettola [1], Laborie [2] has provided a simple but 
expressive formalism for such problems called Resource Temporal Networks (RTNs) 
that we will study in this paper. Briefly, RTNs consist of a Simple Temporal Network 
(STN) as described in [3], constant resource impacts (either production or consumption) 
for events, and piecewise constant resource bounds. Usually, scheduling is performed 
assuming that the problem's characteristics are known in advance, do not change, and 
that the execution of the schedule is deterministic. Often, these assumptions are violated 
in practice. For example, if events do not take place exactly when they are scheduled, it 
may be costly to find a new schedule. Maintaining temporal flexibility during scheduling 
(as was done in knnAEPS) permits the construction of a schedule without determining 
exactly when events take place until execution. Building schedules by ordering events 
rather than assigning event times preserves temporal flexibility. Techniques such as that 
described in [5] make it possible to efficiently update the flexible schedule once the 
precise timing of events are known. 

A simple resource use model assumes all resources impacts occur instantly; this is 
known as the piecewise constant resource impact model. Under this assumption, it is 
straightforward to calculate the resource utilization overtime of events whose execution 
time are fixed. The task becomes more difficult when activities or events are not fixed 
in time, and more difficult still when events can have arbitrary impact on resources. 
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An important use of these techniques is to provide a means for halting the scheduling 
process, either by determining that the resource constraint is satisfied or proving that it 
can’t be satisfied, implying that some scheduling decisions need to be changed. In the 
context of constructive search, early detection of success or failure is often important to 
achieving good search performance. While resource bounds are also important inputs 
to heuristics to drive search in this paper we focus on the ability of resource bounding 
techniques to reduce the cost of constructive search algorithms by means of detecting 
success or failure. 

Two existing techniques address the bounding of resource availability for activities 
or events with a constant resource impact, the Balance Constraint ( BC ) [6] due to 
Laborie, and the Resource Envelope (E t ) [1] due to Muscettola. These techniques are 
described in more detail in Section 2. The techniques are somewhat different, with 
BC featuring an efficient (but loosely) bounding approximation, while E t is somewhat 
more costly but provides a tight bound (in all cases a schedule justifying the bound is 
proved to exist). Somewhat surprisingly, these techniques are not strictly comparable in 
terms of the size of the search trees generated under chronological search. We provide 
examples demonstrating this in Section 3. In Section 4 we describe a generalization of 
these techniques to tightly bound the amount of available resource for a set of partially 
ordered events called the Flow Balance Constraint ( FBC ). FBC is a synthesis of 
the approaches described in [6, 1]. We prove that FBC generates smaller proof trees 
under chronological search. In Section 7 we exploit the results of [7] to calculate FBC 
incrementally, thereby reducing its computational cost. In Section 8 we then generalize 
FBC in order to construct even tighter resource bounds but at increased computational 
cost. Finally, in Section 9 we conclude and describe future work. 


2 Previous Work 


In this paper we will assume that time is represented using integers 1 . We will also 
assume that the resource impact of each event is known when solving begins. We will 
assume each RTN has only one resource, that there is one lower-than constraint and 
one greater-than constraint on the resource, and that the greater-than constraint imposes 
a minimum resource availability of 0. We also assume that the lower bound of the 
scheduling horizon is 0, and that the resource has its maximum availability at time 0. 
The results that follow do not depend on these assumptions. 

The techniques we study are aimed at checking the Necessary Truth Criterion [2, 
8] (NTC), namely, whether there exists a feasible solution of the STN that also satisfies 
the resource constraints. As with all constraints, for a given RTN the NTC can be ei- 
ther proved satisfied, proved unsatisfiablc, or neither. The NTC is proved satisfied if the 
maximum calculated availability of the resource is always an upper bound on the actual 
maximum available resource and is below the upper limit, and the minimum calculated 
availability of the resource is always a lower bound on the actual minimum available 
resource and is above 0. The NTC is proved unsatisfied if the maximum calculated 
availability of the resource is always an upper bound on the actual maximum available 
resource and is below 0, and the minimum calculated availability of the resource is al- 
ways a lower bound on the actual minimum available resource and is above the resource 
upper bound. Otherwise, the satisfiability state of the NTC is undetermined. 

We will use the following notation: Let J\f be the set of all events of an RTN and 
n = [jV’l. Let X 6 c(X) denotes the resource impact of X. If c(X) < 0 X is said 
to be a consumer ; if c(X) > 0 then A" is said to be a producer. As defined in [1], A 
anti-precedes Y if A must occurs at or after Y (i.e. Y < A). A predecessor set of X 
is a set S{X) such that if X e S then every event Y such that Y < A” is also in S . 
A successor set of X is a set T(A’) such that if X G S then every event Y such that 
Y > X is also in S. Let R be an RTN and let A(R) be any procedure for evaluating 


1 This assumption can be relaxed, but leads to resource bounds holding over half-open intervals 
of time. 



the NTC If A(R) — T then the NTC is provably satisfied. If A(R) = F, the NTC 
is provably unsatisfied. If A(R) =? the NTC is neither provably satisfied nor provably 
unsatisfied. 


2.1 The Balance Constraint 

Laborie’s Balance Constraint ( BC ) [6] calculates bounds on the maximum and mini- 
mum amount of a resource immediately prior to and immediately following the execu- 
tion of an event X. The STN is partitioned relative to X into the following sets: Before 
B(X), Before or Equal BS(X), Equal S(X), Equal or After AS(X), After A{X ), 
and Unordered U{X). These sets are then used to calculate bounds on the following 
quantities: £< in (X) for the minimum available resource before X happens, L^^X) 
for the maximum available resource before X happens, L> in (X), for the minimum 
available resource after X has happened, and L^ Lax (X) for the maximum available re- 
source after X has happened. A sample calculation: let P< ax (X) — (B( X) U (V £ 
BS(X) U U(X)\c(V) > 0) Then an upper bound on L< ax ( X ) = T.z&p<^( x) C ( Z )- 
The other bounds can be constructed in a similar manner. The bounds are loose in that 
no schedule consistent with the existing constraints may achieve the bound. The com- 
putational complexity of BC is dominated by the maintenance of arc-consistency of 
the STN. 


2.2 The Resource Envelope 

The Envelope Algorithm E t [l] takes a different approach to bounding resource avail- 
ability. The approach is to find schedules permitted by the STN that maximize or mini- 
mize the available resource at a given time t. The envelope of an RTN is the pair of func- 
tions from time to the maximum (T max (t)) and minimum ( L’minif )) available resource 
defined by the current STN. Events are partitioned into those that are Closed C(t ), 
Pending P(t), or Open 0(t). The maximum available resource at a time t , L max (t), 
is supported by a schedule ensuring that the events in (P m ax(t) C P(t)) U C(t) all 
occur before t . To find the set AP^^ft) C Pit ) that contributes to P~> a:r (£) ? a maxi- 
mum flow problem is constructed using all anti-precedence links derived from the arc- 
consistent STN. The rules for building the flow problem to find L max {t) are ^ follows: 
all events in P(t) are represented by nodes of the flow problem. If X anti-precedes Y 
then the flovc problem contains an arc X — ► Y with infinite capacity. If c(X) > 0 
then the problem contains an arc a — * X with capacity c(X ). If c(X) < 0 then 
the problem contains an arc Y — r with capacity |c(X)j. The flow problem to find 
Lmin(t) is constructed in a similar manner. Examples of the flow problem construction 
are shown in Figure L The maximum flow of this flow network matches all possi- 
ble production with all possible consumption in a manner consistent with the prece- 
dence constraints. Hie predecessor set of those events reachable in the residual flow 
is AP max (t). Define AP£ ax (t) = P(t) - AP max (£). The tightness of the bound is 
guaranteed by proving that adding the constraints {X u b < £}V X € Z\P max (t) and 
{ Y ub > t}VY £ AP^ ax (t) U 0(t) is consistent with the original STN. 


It mms out that AP rnax ( t ) need only be computed at < 2 n distinct times, w ? hich we 
refer to as Instants. The set of these times I — {VXX^ U X u b} is the set of all lower 
and upper bounds of events. Obviously the set of events that could define L max (t) ^ oes 
not change between consecutive instants. We will usually refer to the unique instants. 
We will also refer to e(i) = {i\Xib = i U F|lu6 = i}, those events that define i in /. 
The complexity of the algorithm described in [1] is 0{n * MaxFlow(n,m )), where 
n — |A r j and m is the number of anti-precedence relationships in the arc -consistent 
STN. In [7] this is reduced to O (Max Flow (n : m)) by taking advantage of the order 
in which edges are added to and removed from a single maximum flow problem. The 
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(a) Resource Temporal Network 



<[1,10], -1> <{2,1 1).+1> <C2,11].-1> <3,12].+1> 



Fig. 1. The Flow(s) for E t . 


crucial observation is that when computing E t an event always move from open to 
pending to closed and stays closed. As shown in Figure 1, this guarantees that arcs and 
vertices of the flow problem are removed from the sink r and added to the source cr of 
the flow problem. The insight is that the flow problem need not be solved anew, but the 
previous flow can be reused, thereby saving a factor of n. As an aside, a more general 
incremental maximum flow algorithm is given in [9], but it doesn’t take advantage of 
the order of flow arc additions and deletions, and is not as efficient at calculating the 
envelope. 

3 Examples of Non-Domination 

Muscettola previously demonstrated that E t can prove the NTC is satisfied in cases 
where BC cannot. In this section we provide an example where E t fails to prove the 
NTC is violated in cases where BC can. The somewhat surprising consequence of this 
is that neither algorithm ’’dominates” the other when used solely to halt search. 

3.1 BC Doesn’t Dominate E t 

Muscettola [1] describes an RTN for which BC does not prove the NTC is satisfied. 
This example is modified and reproduced in Figure 2. Initially the resource has 2 units 
available. BC will reason as follows: U(A) = {W.X, Y, Z} and A(A) = { B , C, D}. 
Then L^ in (A) = -1; it schedules W and Y before A, and schedules X and Z after 
A. Clearly, this schedule is not consistent with the original STN, so the bound calcu- 
lated by BC is ’’loose”. Thus, BC incorrectly believes that more decisions are needed. 
The E t algorithm is able to prove that the NTC is satisfied in this case by determin- 
ing L m i n (t) > 0 over the scheduling horizon. As a case in point, consider L rn i n ( 11). 
The pending set P(ll) = {B, C, D, X, Y, Z}. We see that it is possible to postpone B 
and X and achieve L min (ll) = 0. It is also possible to schedule C and Y , the other 
consumers; however, scheduling C, for example, forces the scheduling of the producer 
B prior to C because of temporal constraints. This accounts for one of the consump- 
tions, and makes it impossible to improve upon L m i n (ll) — 0- If we consider the 
flow problem in Figure 1 (b), we see that this is exactly what is calculated. In addition. 
Lminit) < 2 over the horizon, showing that the NTC is provably satisfied. 



(a) Resource Temporal Network 
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(b) Envelope and Balance Constraint lower bounds 


Fig. 2. An RTN for which BC fails to detect that the problem is solved. 


3.2 E t Doesn’t Dominate BC 

Figure 3 describes an RTN for which E t cannot show that the NTC is provably not 
satisfied. Again, initially the resource has 2 units available. In this case, the BC will 
find L< ax (A) = L< ax (B ) = L< ax (C) = 1, but L< ax (D) = -5, since B(D) = 
{A,B,C}. This proves that, no matter what additional constraints are imposed, the 
NTC will be violated. By contrast, E t cannot conclude that the NTC is violated. We see 
that B and C can be postponed until time 10, and A can be scheduled early, leading to 
L-maxQ) — 3- At time 2, we could schedule D, but that would require scheduling ev- 
erything before 2 with a resource availability of 0 being the result; at time 2 we can still 
find a schedule in which A and nothing else so L rnax (2) = 3. At time 10, however, both 
consumption events must have taken place; in order to achieve the maximum resource 
available we can schedule D, leading to L max {ll) — 0. The calculation of 
(not shown) provides no assistance in proving the NTC is violated. Not only does E t 
fail to show that the NTC is provably violated, there are non-trivial ordering decisions 
that can be made; in the worst case, schedulers could spend a considerable amount of 
time continuing the search from states like the one shown in Figure 3. 


4 A Tighter Balance Constraint 

4.1 Setup 

In order to achieve tighter bounds than either E t or BC , w^e adopt a synthesis of 
both strategies. We use the following partition of the STN relative to X : C( X) is 
the set of events that must have been scheduled before X, 0(X) is the set of events 
that must be scheduled after X, E( X) is the set of events that must occur at the 
same time as X, and P(X) is the remainder of the events. Like Laborie, we then find 
L max( x ): L min( x )’ L max( x ); and L min( x )• Like Muscettola, we build a maximum 
flow problem wdiose residual graph is used to construct the supporting sets P< ax (A"), 
PP n {X). P> ax (Xl and PC n (X). The rules for constructing the flow problem are 
identical to those described previously; only the set of events considered in the flow 
problem is different. The resulting set P r ^ ax (X) defines an STN constructed by adding 
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Fig. 3. An RTN for which E t fails to detect that the problem is unsolvable. 


the following constraints: W e ^ai(^){^ < X} and VV € P^ x (X){ X < V}. 
Note that we also need find AP< ax (X) C P(X) and so P^ ax ( X) C P{X) U C( X). 
The other set containment properties similar. We refer to the resulting bounds as the 
Flow Balance Constraint (FBC), since it combines the features of the envelope with 
the flow-based approach to guarantee tightness. 

One might think from this analysis that the number of flow problems to solve is 4 n, 
i.e. one flow problem for each of L^ in (X),L^ in (X),L< ax (X),L% lax (X). Some- 
what surprisingly, this is not necessary; only 2 n flow calculations are needed. Consider 
P(X). By definition, X has no anti-precedence arcs to any event in P{X); in the lan- 
guage of [1] it is flow-isolated. X doesn’t contribute to L^ lin (X) 7 L^ iax (X) and does 
contribute to L^ lin (X),L^ iax (X); once the flow problems are solved, we simply add 
E(Xy s resource impact. The sets we find from solving the flows are thus referred to as 
AP^JX) and AP^ in (X) 


FBCOl) 

Make the STN S arc-consistent. 

Infer all anti -precedence and precedence links, (transitive closure) 
Collect all sets P(X), C(X). 
for each event X 

Build flow problems for upper and lower bounds from P(X) 
Bound(X) 

end for 

end 

Bound(X) 

Find AP^ax (X) C P{X) 

Lmax(X) — 5Z Ve / i p« iax (x)uC(X)uE(X) C O'O 

Lmax(X) = ^2 VeA pC ia ^^ X )UC(X) 

Find AP^iniX) C P{X) 

^min(X) = y^veAP^ n {X)UC{X)uE(X) C W) 

= Ylv€AP^. n (X)L\C(X) ) 

end 


Fig. 4. A sketch of the Flow Balance Constraint. 



The algorithm for FBC is described in Figure 4. Since we must derive the transitive 
closures of the anti-precedence and strict precedence constraints, we can collect the sets 
P{X), E{X), C{X) in time proportional to m the number of induced anti-precedences 
in the graph. Because we find AP^ ax (X) and AP^ in (X) and using maximum flow, 
this is arguably the most expensive part of the algorithm; in the next section we discuss 
complexity in more detail. 

We now proceed to prove that the resulting bounds on quantities like L^ iax (E) are 
tight To do so, we first show that there is at least one schedule justifying the bound, 
and then show that there is no better bound than that found using the flow problem. 

Theorem 1. Let R be an RTN and X be an event in R. Suppose AP^ iax (X) ^ 0. 
Let R' be the STN formed by adding the following constraints to R: {V < X} for all 
X e AP^ ax (X) and {X > V} for all X e AP£ ax (X). Then R' has at least one 
temporally consistent solution . 

Proof Since V € P{X) the imposition of a single constraint alone doesn't make the 
STN inconsistent. Imposing a constraint V < X only decreases V and increases 
Xu,. Since AP^ ax (X) is a predecessor set, all {V < X} can be imposed simulta- 
neously without impacting consistency. Imposing a constraint X < V only decreases 
X ub and increases V Lb . Since AP£ ax (X) is a successor set, all {X > V} can also 
be imposed simultaneously without impacting consistency. Finally, X is the only event 
whose bounds are acted on by both classes of constraint. Consider two events A,B. 
Since A € AP^ ax (X) and B e AP£ ax (X) we know it can't be the case that B < A . 
But then either A < B or A and B can be ordered in any way, and we already know X 
can be ordered any way with respect to A or B. Thus, A < X < B is possible and no 
such ordering prevents other linearizations with respect to X. □ 

Theorem 2, Let R be an RTN and X be an event in R. Suppose AP^ ax (X) 0. 
Then HveAP^jx)uC(X)uE(X) c ( v ) is the maximum possible value ofL^JX). If 
^ZxeE(X) C ( X ) < 0 l ^ en TIv^ap* ( x)uc(x) C (Y) ^ ie maximum possible value 
for L< ax (X), otherwise L^ ax (X) is the maximum possible value for L^ ax (X). 

Proof Since we construct the flow problems in exactly the same way as is done in [1], 
we state this as a corollary of Theorem 1 of [1]. □ 

The proofs for the tightness of the bounds on L^^f X). LF n (X) and L^ in (X), 
are similar. 


4,2 Relating the Envelopes 

In this section, we formally establish the relationship between the intervals over which 
the FBC and E t bounds hold. 

Since the envelope is the resource bound as a function of time, we should be able to 
determine the interval of time over which the event based bounds hold; then the max- 
imum (minimum) envelope at t is the upper (lower) bound of all of the event based 
bounds that hold at t. Intuitively, L^ iax (X) is the maximum availability of the resource 
for an interval of time immediately prior to the time X happens. What is the interval 
of time over which this value holds? Suppose we find a set P< ax (X) that supports 
L^ iax (X). Then as previously stated, this bound assumes that we impose the follow- 
ing constraints: VV e P< ax (X){V < X} and YV € P^ x (X){X < V} to get 
a new RTN TV. Note that the new pending set P f {X ) = 0 and the new' closed set 
C ( ( X) = P^ ax (X). Note also that L^ ax (X) is the resource availability 7 no earlier 
than V* b , = max vec'(X) - Due to the new constraints X may have a new up- 
per bound X ub > , which defines the latest time that L^ ax (X) holds. Finally, note that 
L max ( X ) = L max' (t) over the interval [V* b , , X ub A y 

We now can demonstrate the precise relationship between L max (t) and L^ ax ( X). 


Theorem 3. L^ ax (X) < L max (t)for V* b , < t < X ub >. 


Proof. Since we add new constraints, 'itL max >{t) < L max (t). Since L^ ax (X) = 
Lmax'{t ) we’re done. □ 

Now suppose we find a set P^ ax (X) that supports L^ ax (X). Then as previously 
stated, this bound assumes that we impose the following constraints: W e P> ax (X){V < 
X} and iV 6 P>£ x (X){X < V'} to get a new RTN 7?/. Note that the new pending 
set P'{X) = 0 and the new closed set 0'(X) = P> ax (X). Note also that L^ ax (X) 
is the resource availability no later than Wf b , - minvi/ e c'(X) W ib -. Due to the new 
constraints X may have a new lower bound Xw , which defines the latest time that 
L max( x ) holds. Finally, note that L^ ax (X) = L ma x' (t) over the interval [X lb > , Wf b ,\. 

An example is shown in Figure 5. 

Theorem 4. L^ ax (X) < L max {t) for X w < t < Wf b ,. 

Proof. Identical to the previous proof. □ 


(a) Resource Temporal Network decomposed relative to event X. 



(b) The new Resource Temporal Network. The bound v* is the latest time of any 
event in P^^X). L max > (A) holds over the interval [V* ufc) ,, X ub ,]. 


Fig. 5. Deriving the time intervals over which resource availability bounds hold. 


We show an example of the relationship in Figure 6. The RTN is the same as that 
in Figure 3. This example also shows why BC can answer positively in cases where 
E t cannot. In this case, BC finds a tight bound on L^ iax (D) that proves that the NTC 
cannot be satisfied. Unfortunately, E t must calculate the maximum over all of the event- 
based bounds that hold at a time t , and cannot do better than return ? for this RTN. 

5 Dominance 

In this section we will describe our dominance criteria and show that FBC dominates 
E t and BC: we do not formally show that BC and E t do not dominate each other, re- 
lying on the intuition of Section 3 to adequately demonstrate this. As previously stated, 
we assume that FBC. E t and BC are primarily used to detect whether a scheduling 
algorithm can halt or must continue. In order to motivate our definition of dominance, 
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Fig. 6. An RTN for which E t fails to detect that the problem is unsolvable. 


suppose that some resource bounding procedure A is used in a chronological search 
framework with a fixed variable and value ordering heuristic Then we would like A 
to domiante B if A leads to smaller search trees than using B . Now suppose A and B 
are used in a stochastic sampling approach where we sample RTNs that are temporally 
consistent with the original problem. Then we would like A to domiante B if A leads 
to less sampling than B. We use the following definition of dominance. 

Definition 1. Let R be an RTN. Let T(R) be the set of all temporally consisten t RTNs 
that can be formed from R by adding temporal constraints to R. Let A B : R —+ 
{T 1 , F, ?} Let U a (R ) = 5 € F(i?)|,4(S) =?. Then A dominates B on R ifU A {R)) C 
Ujs{R)). *4 dominates B if3R such that A dominates B on R and there exists no S 
such that B dominates A on S. We write A ~<r B or A A B as appropriate. 

Theorem 5. FBC A BC. 

Proof Theorems 1 and 2 show that the flow construction guarantees the tightest possi- 
ble bounds on L< ax (X). L> ax (X),L< m (X) and £> in (X) for any RTN. Thus there 
can be no S such that BC A 5 FBC. The example shown in Figure 2 shows at least 
one RTN R for which FBC A# BC. This completes the proof. □ 

Define E<(t) - X\L< ax (X) holds at t and E>(t) = X\L> ax (X) holds at t. 

Theorem 6. L max (t ) = max(max xsE < (t ) L< ax (X),max XS E>(t) L-> ax (X)) and 
Lminit) = min(min X €ij<(t) L< in (X),mm X€B > (t) L^ in (X)). 

Proof Corollary’ of Theorems 3 and 4 □ 

Theorem 7. FBC A E t 

Proof Theorem 6 shows that there is no RTN S for which E t As FBC. The example 
in Figure 3 shows that there is at least one RTN R for which FBC A# E t . This 
completes the proof. □ 

6 Complexity 

A naive algorithm for calculating FBC builds a new flow network for each event 
X € A r . An equally naive analysis of the complexity of this algorithm is 0(n * 
M axFlowfn. m)), where rn is the number of anti-precedence relationships in the RTN. 
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However, this is unsatisfactory for a number of reasons. Consider the RTN in Figure 
2. The flow networks required to calculate FBC are identical for A, B^C.D because 
P(A) = P(B) — P{C) = P(D ) . Additionally, the flow networks include only a small 
fraction of the n nodes in the RTN; strict precedences and equalities among events will 
generally reduce the size of P{X) 9 leading to smaller flow problems. These observa- 
tions suggest it should be possible to analyze the structure of the anti-precedence graph 
and reduce the number of flow problems to solve, with a resulting tighter bound on the 
complexity of calculating FBC. 


In this section we provide a lower bound on the complexity of the naive approach 
to calculating FBC. We do so by constructing an RTN such that the flow problem to 
solve for each event is both non- trivial and distinct. The RTN is a ’’square” graph with 
y/n events per side. We index events by row and column in the square. The RTN has 
the following strict precedences: (i, j) < (i 4- 1 ,j) 9 an< ^ (hj) < (hJ + 1) (obviously 
omitting those links for which the indices are outside the bounds [0, y/n\. By construe - 


THfi ;\\ — ff~ --M- 

MUjj. u °yJJJ — VA^A^F 


> ... ^. *\ ■ , (f M... 

* y ^ J) ^ u a » v )i u 


■ * » • N rfiT , . _ T 1 _ -C a! . _ 

l / \ U > j). lUUij, <111 Ui U1C uuw 


graphs are distinct. Every such set has a non-trivial part of the flow graph. Notice that 
we can assign c(X) arbitrarily to the events of the RTN, as long as they are all non-zero 
and there is a mix of consumers and producers. This RTN is shown in Figure 7. 




(b) The pending set for node (3,3). 



Fig. 7. A perverse RTN providing a worst-case lower bound for calculating FBC. 


We now proceed to construct a lower bound on the complexity of the naive approach 
for calculating FBC on this RTN. By construction we have guaranteed that no ’’quick 
fixes” can be used to decrease the complexity. The larger of the two induced flow prob- 
lem forevent (i,j) contains max((i - 1)(\/^ ~ j — 1). (j - 1 )(y/n — i — 1)) events, 
and at least this many flow arcs (we could do an exact count but it isn’t necessary since 
we’re providing a lower bound.) Let us now assume we are using a FIFO preflow -push 
algorithm to solve each flow problem [10]; this ignores any efficiency gained from an- 
alyzing the pushable flow, but is also suitable for our purposes. If v is the number of 
nodes in the flow problem, there are at least v edges, and so the complexity is i?(?; 2 °). 
Using these assumptions the total complexity of solving all the flow problems is 



y/n \/n 

y^y^ max((i - l)(Vn - j - 1), (j - l)(y/n ~ i - l)) 2 ' 5 

J — 1 i— 1 

First, we simplify the sum to get a lower bound: 

\/n— 1 t/ti— 1 

< E E Wv^-j )) 2 - 5 

jf=0 z— 0 

We next approximate the sum with the integral (which, while bounding above, is 
close enough for our purposes): 

/ vn-l / rVn-1 _ \ 

( J_ ~ j)) 2 5 di J dj 

The first integral with respect to z is trivial. Substituting for y/n ~ z, we see that 
1 (v^ ~ — (V™ “ z) 2 '° j , so ultimately we get 

„ / (yn- 1) 35 \ / y/n 3 ' 5 1 \ 

\ 3.5 ) 3.5 3.5 J 

Collecting the high order positive powers of y/n we see that the naive algorithm has 
a lower bound i2(n 3 °). Thus, for preflow-push flow algorithms using FIFO queues, the 
naive algorithm for FBC is Q(n * MaxFlow(n , m)). 

7 Incrementally Calculating FBC 

As stated previously, the incremental technique described in [7] shaves a factor of n off 
the cost of calculating F t . This provides some hope that we can find a way to eliminate 
the factor of n ’’extra’* cost for calculating FBC . A naive approach to employing the 
results in [7] for calculating FBC is inadequate to eliminate a factor of n from the 
complexity. The crucial element of the complexity’ analysis in [7] requires that an event 
always move from open to pending to closed. We show that naively applying the incre- 
mental algorithm may result in a non-trivial number of events moving from closed to 
pending, thereby defeating the cost- savings measures. Consider the RTN described in 
Figure 7. The longest chain of events for w 7 hich the incremental algorithm can be used 
to calculate FBC and reuse the previous flow 7 is of length 2 y/n (either the diagonal or 
two edges). Each of these induces a total complexity 7 of 0(M ax Flow (n. m)) because, 
eventually, all events and arcs in the RTN will be considered in some now problem. 
There are 0{y/n) such chains; once one chain is done, all the information about the 
previous flow must be reset, because at least one event that w 7 as previously closed be- 
comes pending or open. The total complexity is therefore 0(y/nMaxFlow(n, m)). 
At worst, no ordering of the events can be found to eliminate this; at best non-trivial 
analysis of the precedence graph is necessary to find such an ordering. 

A solution to this problem involves using the topological structure of the precedence 
graph to order the calls to Bound(X) and to cache the information from each flow 
problem that is solved. As long as E C(X ,:) no event in C(Xi — 1) ever enters 
P(X t ) and no event in P(Xi-i) ever enters 0(X l ): this ensures that the incremental 
approach of [7] can be used. On the other hand, whenever X t does not strictly follow 
i, w T e have cached the necessary' data from the flow problem for some event Y 
that strictly preceded X^ again ensuring that the incremental approach of [7] can be 
used. The node order requires beginning with a node with no predecessors. Implicitly 
the node order traverses the precedence graph using a topological sort. It’s trivial to 
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initialize the set of events with no precedences, and also trivial to mark each visited 
node to ensure a topological order is obeyed. 

For the analysis of the algorithm, we simply fix an arbitrary topological order and 
concentrate on describing the storage of the flow r information and the effective size of 
the largest flow problem that could be solved. We call the resulting algorithm FBC — 
DFS , and it is described in Figure 8 below. 


FBC-DFS(T^Xi) 

Make the STN S arc-consistent. 

Infer all anti-precedence and precedence links, (transitive closure) 
Collect all sets P(X), C(X). 

Collect set y of events with no predecessors 
F-max “ Fmin — 0 


for each event Xi 


if ^Xi^eC(Xi)) 


V \ 


Pop Fmax and Fmin to Y 

end if 

Build flow problems for upper and lower bounds from P(X),top(F m ax) and top(Fmin) 
Bound(Xi) 

Push the flow solutions and events they correspond to onto F maI and Fmin 

end for 


end 


Fig. 8. A sketch of the Flow Balance Constraint implemented by using a Depth-First search 
through the precedence graph. 


Theorem 8. Let 7Z be an RTN with m induced anti- precedence constraints in its STN 
S. Algorithm FBC — DFS takes 0(MaxFlow(n , m)) time and 0(nm) space. 

Proof Ensuring events are searched in topological order imposes no significant over- 
head, nor does finding the correct event to back up to. If Xi strictly precedes X in the 
order that Bound(A'j) is invoked, then the incremental complexity argument of [7] ap- 
plies. Essentially, we must (eventually) construct the flow problem for the entire STN. 

At worst, each time we call Bound(TQ) we must store the resulting flow network. 
Because we are searching the graph depth-first, at worst we do so O(n) times. If not, 
we must pop the stacks Fmax and Fmin each time we ’’back up”. At worst we back up 
0(n) times and store data on each of m flow edges. Only a constant amount of data is 
required for each edge. This cost is O(mn), which is dominated by the complexity of 
the flow (for FIFO prefiow-push, anyway, because yTn < n). □ 

Note that weaker conditions on the events may lead to improved computational 
performance, but the conditions we have imposed are sufficient to produce an algorithm 
with the desired run-time complexity. 


8 Higher Order Balance Constraint 

The quantities L< ax {X), L> ax (X), L< in {X) and L> in (X) can be thought of as first 
order checks on resource availability, in that they calculate resource bounds before and 
after one event. In this section we generalize these techniques in order to calculate 
higher-order resource availability checks after k events. To see why this is valuable, 
consider Figure 9. We see that no first order bound calculated by FBC proves that 
this RTN is impossible to solve. However, notice that we can show that the maximum 
available resource after both A and B have occured is —1. This corresponds to one 
of two schedules: it is possible to either schedule D < B or C < .4, but not both 
simultaneously. Thus, without further search, we can prove that the NTC is violated. 
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(a) Resource Temporal Network 1 st and 2d order Flow Balance Constraint 

upper bounds. The second order bound 
L max > (A A B) proves the NTC is violated. 


Fig. 9. An RTN for which FBC fails to disprove . 


This example shows that it may be valuable to perform 2^-order checks on resource 
availability by determining resource availability immediately before and after sets of 2 
events. In order to do this for L^ ax (X A F) we must account for the following possi- 
bilities: neither X nor F have occurred, X has occurred but Y has not, and vice versa. 
First, we solve the flow problem over the events of P(X) n P(F) to find the maximum 
availability strictly before both X and Y . We call this set of events P< ax (-'(x a y)). 
The second requires adding the constraint A' < Y and then solving the flow prob- 
lem defined by P(F) and find the maximum availability strictly before Y. Call this set 
P^a X ( X < F). The last requires adding the constraint Y < X and then solving the 
flow problem defined by P( X) and find the maximum availability strictly before X . 
Call this set of events F< ax (Y < A). We define L^ ax (X A Y ) as 

max( E c(F), 

V€P< aI (^(Xf.Y)) 

E c ( v t 

V€P< ax (X<Y) 

E < v ) 

vep< ax ( x<y) 

(1) 

For L^ ax (X AY) do the following. First, we solve the flow problem over the events 
of P(A~) U P(F). Call the resulting set P> ax ((X * F)). Now we must also include 
events in E( X) U P(F), but some of these may already be in P(A) U P(F), so these 
events must be discarded. Define R( X V Y) ~ (E( X) U E{Y)) P (P(A) U P(Y)) U 

P max(( X * Y ))- We A 05 define L m a x(( X V Y)) = Ev€i?(XvK) C ( V )- This is the 
amount after the schedule assuming both X and Y have occurred. The lower bounds 
are calculated in a similar manner. 

To see how this works in Figure 9, note P(A) PP(P) = 0. The maximum available 
under these circumstances is achieved by postponing all the events, w^hich is 2. If we 
impose .4 < B, P(B) — D ; the best schedule here is A, D for an availability of 1. 
The same applies if B < A. Thus, L^ a JX A Y) = 1. To calculate L> cax (A A B) 
we observe that the first incremental calculation leads to 2 — 4 = —2 since neither 
A nor B are assumed to have occured. If we impose A < B, the max involves A ; 
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but we must assume B has happened as well to achieve the maximum, which leads to 
1 — 2 = —1. The same applies for impose B < A. The result leads to the schedule 
defined by ensuring all of the events are before or equal to A and B. We can define the 
interval over which the bounds hold as we have previously. 

The total complexity of the resulting naive algorithm for calculating Second order 
Flow Balance Constraint ( FBC 2 ) is Q ((n 2 ) Max Flow (n,m)). The n 2 term comes 
from the fact that n(n - 1) pairs of bounds must be calculated; the complexity bounds 
obscure the fact that 3 flow problems must be solved per bound. 

Note, however, that n 2 is a very crude estimate of the total number of bounds to 
compute. If A strictly precedes B then P^ ax (A A B) = P< ax (B). Thus, the induced 
precedences vastly reduce the number of bounds to calculate. Additionally, the sizes of 
the flow problems will generally be larger as the number of events involved climbs. This 
is because P(X A Y) = P(X ) U P(Y). These factors make a more precise complexity 
analysis difficult. Finally, it is likely that the incremental flow algorithm described in 
[7] can be used to further reduce the complexity. It is sufficient for our purposes to 
demonstrate that even tighter inferred constraints can be calculated in time polynomial 
in the number of events considered. 

We can further generalize this to sets of k events in the same manner. The complex- 
ity of the naive algorithm for calculating FBC k is Q (2 fc ) MaxFlow(n. m)^ . 


This is because there are 


bounds to calculate, and each bound requires solving 2 k 


flow problems (as well as arc-consistency enforcement steps). 


9 Conclusions and Future Work 

In this paper we have shown, contrary to expectations, that BC and E t are not strictly 
comparable in terms of their power to halt search over partial ordered schedules for 
RTNs. We have also shown how to exploit the features of BC and E t to construct 
FBC , a tighter bound on the availability of resources for RTNs than either of the previ- 
ous approaches. The resulting bound can be computed by the algorithm FBC — DFS 
in 0{MaxFlow(n : m) time, but n * O (Max Flow (n,m) space. When used in iden- 
tical search algorithms with identical static variable and value orders, FBC will gen- 
erate search trees with less than or equal nodes than either BC or Ft- The technique 
generalize for calculating FBC leads to even tighter bounds, but at sharply increased 
computational cost. 

While we have proven theoretical dominance of FBC, an empirical study will be 
necessary to determine whether it is worth the overhead. An empirical study will also 
have the added benefit of shedding more light on the relative value of E t and BC for 
speeding up the solution of scheduling problems. There are numerous possibilities for 
speeding up the solution of the flow problems necessary to calculate the bounds. Fur- 
thermore, despite their apparent complexity, higher order variants of FBC may also 
prove worthwhile. The DFS-FBC algorithm could benefit from judicious node order- 
ings. One possible reason to order the flow problems is to reduce total flow problem 
solving costs, another (possibly conflicting) goal is to reduce storage costs. Similarly, 
the higher-level consistency algorithms can be crafted to maximally exploit the incre- 
mental envelopes approach as well as judicious ordering of the sets of k events that 
must be contemplated to enforce the bounds. 

Changes to dominance criteria that we use are worth contemplating. On the one 
hand, requiring that A dominates B on R if Ua(R)) C Ub(R)) is rather strong, and 
could be weakened, say, to \U^(R))\ < |f7g(i?))j. On the other hand, requiring that A 
dominates B if 3 R such that A dominates B on R and there exists no S such that B 
dominates A on S is somewhat weak, and perhaps could be strengthened. 

A second, equally important empirical study will be necessary to shed light on how 
to integrate heuristics that make use of the various bounds. Laborie [6] has built numer- 



ous such heuristics for BC , providing a good starting point. However, such heuristics 
likely have complex interactions with the pruning power of the envelopes. It will likely 
be necessary to trade off between the pruning power and heuristic predictiveness of the 
resource bounds to craft the best scheduling algorithm. 

For simplicity, we concentrated here on analyzing the complexity of FBC using 
Prefiow-Push with FIFO queues. Generalizing the complexity analysis to other flow 
algorithms may be worthwhile, but is complicated by several factors. The pushable flow 
cannot be easily analyzed, but it is tempting to try assessing the relationship between 
the number and size of the flow problems to be solved during calculation of FBC by 
looking at the precedence and anti-precedence graphs. Appealing to graph theory may 
lead to both improved complexity analysis and, possibly, better algorithms as w'ell. 
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